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Abstract 

We explain how static multi-vortex solutions arise in non-linear field 
theories, by taking the non-linear Schrodinger equation coupled to Chern- 
Simons field (Jackiw-Pi model) and a fermion Chern-Simons theory as 
simple examples. We then construct a fermion Maxwell-Chern-Simons 
theory which has consistent static field equations, and show that it has 
the same vortex solutions as the Jackiw-Pi model, but gives rise to quite 
different vortex dynamics. 



*To appear in the proceedings of "International Workshop on Mathematical and Physical 
Aspacts of Nonlinear Field Theories" , 1998. 
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I. Introduction 



Magnetic vortices are found to appear in high temperature superconductors as 
quantized vortex lines in the range of magnetic field between 10 -2 tesla and 
10 2 tesla [0. It is known that the behavior of vortices dominates many phys- 
ical properties of high temperature superconductors in this range of magnetic 
field. These vortices show very interesting dynamical properties:in some range 
of magnetic field and temperature they form vortex lattice or behave as liquid. 
This phenomenon poses an important theoretical challenge in understanding 
the dynamics of magnetic vortices. 

One way to understand the dynamics of vortices is to consider non-linear 
field theories that possess classical vortex solutions and to study the dynamics of 
these solutions in the field theoretic framework. There exist many field theoretic 
models that support the static vortex solutions. Relativistic models of vortices 
include the Abelian Higgs model of Niclson and Olesen the scalar Chern- 
Simons theory of Hong, Kim and Pac, and Jackiw and Weinberg ||, and the 
fermion Chern-Simons thoery of Lie and Bhaduri Q. There also exist many 
non-relativistic models including the Ginzburg-Landau model || which is the 
non-relativistic limit of the Abelian Higgs model, Jackiw-Pi model Q which is 
the non-relativistic limit of the scalar Chern-Simons theory, and non-relativistic 
spinor Chern-Simons theories |Q. 

The most of these models are not soluble in closed form except for the Jackiw- 
Pi model and some of its generalizations. For the study of vortex dynamics in 
the field theoretic framework, however, it would be more convenient to have 
completely soluble models. We therefore consider the field theoretic models 
which are completely soluble. 

In the next section we give a brief review of the Jackiw-Pi model to show 
how to find the static self-dual solutions. In section III, the simplest fermion 
Chern-Simons theory which supports the static vortex solutions as solutions 
of the Liouville equations is presented. In section IV, we present a fermion 
Maxwell-Chern-Simons theory which supports static vortex solutions. We show 
that although the static solutions are the same as those of Jackiw-Pi model 
and the simplest fermion Chern-Simons theory, this model gives rise to quite 
different moduli space dynamics from those simple models. We conclude with 
some discussions in the last section. 

II. Non-linear Schrodinger Equation coupled to 
Chern-Simons Gauge Field 

As a simplest model field theory that possesses vortex solutions we consider 
the non-linear Schrodinger field theory coupled to a Chern-Simons gauge field 
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described by the Lagrangian Q, 



C = -^F^ + i<t>*D <f> - \\D<P\ 2 + ^(<^) 2 , (1) 



where 



Ffj,u — A v — d v A^ 
A^ = (A°-A) 

D a = d + iA° 1 D = V-iA. 

One way to see the existence of static solutions is to consider the energy func- 
tional and to find the field configurations that minimize the energy functional 
[gj. To this end it is convenient to write the Lagrangian ([!]) as 

C = ^,,.l'.i< + up* - - kF 12 ) - ±\D<f>f + ^(0*0) 2 , (2) 

where indices i and j run for the spatial components 1 and 2. The momentum 
conjugates to the field variables A 1 and <f) are defined by 

5C k ■ 

7T,' — — = —CjjA J 

5A* 2 13 

P=g = ^, (3) 

0(f) 

respectively. The Lagrangian (||) can then be written as 

C^^A 1 +pj>-H-A°{4>*4>-KF 12 ), (4) 

where H, is the Hamiltonian density and the last term represents the Gauss' law 
constraint, 

B = Vxl= -F u = --<t>*(f>. (5) 

K 

The Gauss' law constraint (Q) shows that the magnetic field is proportional to 
the charge density of the scalar field, 

P = 0*0- (6) 

From Eqs. (||) and (||) one finds that the Hamiltonian density of the system 
is given by 

n= l -\D4>\ 2 -^{4>^) 2 . (7) 

Using the identity, 

\D(j>\ 2 = \{D 1 - iD 2 )(/>\ 2 — (Bp + V x J), (8) 
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where J = Im(f>*D(j) is the current density, one can write the Hamiltonian 
density as 

U = \\{D 1 ~iD 2 ) ( j ) \ 2 -^ x J. (9) 
The energy of the system then becomes 

E = J d 2 xH = d 2 x[\(D 1 - iD 2 )(f>\ 2 - V x J] . (10) 

If the fields are well-behavecd at infinity, then the last term does not contribute 
to the energy since it is a surface term. Then the energy of the system is positive 
definite: 

E = ^Jd 2 x\(D 1 -iD 2 )0\ 2 >O, (11) 

where the equality is called the Bogomol'nyi bound. Note that the coupling 
constant of the quartic interaction of the scalar field is so chosen that the energy 
functional is positive definite. 

Eq. (11) shows that the minimum energy configurations of the system are 
determined by the static first-order differential equation, 

(£>i - iD 2 )(t> = 0, (12) 

which is called the self-dual equation. 

To find the static soliton solutions one has to solve the self-dual equation 
( p"2| ) together with the Gauss' law constraint (0). To do this we write the scalar 
field as 

c/> = ph i ». (13) 
Substituting @ into (|l2|) one finds 

A_ = -i<f>- l d_(j) = -^d-hip + d-uj, (14) 

or 

1= Vuj - -V x lnp, (15) 

where A± = A\ ± iA 2 and d± = d\ ± id 2 . This shows that the gauge field A is 
completely determined by the scalar field. From (Q) we obtain 

d-A + - d+A- = 2iV x A = 2iB = iV 2 In p, (16) 

which, upon substituting into the Gauss' law constraint, reduces to the Liouville 
equation, 

2 

V 2 lnp= p, k>0, (17) 

K 

which is completely integrable. 
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The simplest solution of the Liouville equation (|17|) is the spherically sym- 
metric solution, 

^ = —[(7) +(-) ] . (18) 

where ro and are the constants representing the scale and the flux number of 
the soliton, respectively. To find the restriction on the number N, we observe 
that the regularity of the solution at the origin and at infinity requires N > 1. 
And the single- valuedness of the scalar field <j> requires N to be an integer. Note 
that the regularity of the gauge field at the origin, 

Mr) - \ei~(N - 1), (19) 

determines the function u>: 

u=(l-N)6. (20) 

We thus find the static solution, 

^ = ^[(^)" + (^)"]->-»", (21) 



with A determined by Eq.(|15|) 

From the explicit solution ( ^l| ) one finds the magnetic flux of the soliton 
solution, 

1 



$ = / d 2 xB = --J d X P = -2tt(2A0, (22) 
the electric charge, 

Q = (23) 

and the angular momentum, 

J = 2ttk(2A0. (24) 

This shows that the solution (Q) and ( pl| ) represents the charged vortex solution 
with quantized magnetic flux. 

The solution (|l8|) represents the case of N solitons superimposed at the 
origin. The general solution of the Liouville equation is also known: 



(i + l/WI ) 

where z = re t9 and f(z) is an arbitrary function of z such that p is non-singular. 
f{z) for the most general solution can be written as 

N 

i=X 1 
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where Ci and Z{ are arbitrary constants. This solution contains AN arbitrary 
parameters which represent 2N position parameters, N scale parameters and 
N phase variables of the solitons. In fact, index theorem confirms that (4iV) is 
the maximal number of parameters contained in the general solution || . 

The Jackiw-Pi model described above is invariant under the Galilean trans- 
formation, and is the non-relativistic limit of the relativistic scalar Chern-Simons 
theory § , 

Crei = -e^F^ + D^D^* - 4l0| 2 (H 2 - (27) 

4 ' K 

This theory possesses two types of soliton solutions: the topological solitons in 
the symmetry broken sector, and the non-topological solitons in the symmetric 
sector. The magnetic flux of the topological solitons is quantized as 

$ = 2tt(7V-1), (28) 

while the magnetic flux of the non-topological solitons is not quantized: 

$ = 2n(N + a), (29) 

where a is an arbitrary parameter such that a > N. 

The non-relativistic limit of the model ( p7[ ) in the symmetric sector reduces 
to the Jackiw-Pi model. Why then the non-relativistic solitons have quantized 
flux while magnetic flux of the relativistic case is not quantized? The answer 
lies in the fact that the Jackiw-Pi model has inversion symmetry which is not 
respected by the relativistic model || . To see this note that the charge density 
satisfies the equation, 

V 2 lnp= ( ^f^ 1- £)' relativistic case (3Q) 
\ —\p, non-relativistic case. 

One can easily show that the Liouville equation is invariant under the inversion 
transformation: 

r -> -, 6 -> 6 

r 

p(r) - = r 4 pM, (31) 

while the relativistic equation is not. The charge density for both models be- 
haves as 

, v J r^ N -V, as r^O 
Pto — | r - a(a+ i) j as (32) 

and thus the magnetic flux is given by 
- $ = [ A-dl 



G 



= / (-fw + iv x In p) ■ dl (33) 
= 2n(n - 1) + 2ir(a + 1). 

Due to the inversion symmetry (pl[) of the Jackiw-Pi model, the behavior of 
p(f) at the origin and at infinity (|32[) must be related. This give rises to the 
relation a — N, and explains the quantization of flux numbers in the case of 
the non-relativistic theory. 



III. Vortex Solutions in a Fermionic Chern-Simons 
Theory 

Vortex phenomena in realistic systems are known to be dominated by the elec- 
tronic structure of the systems. Thus in attempting to understand the vortex 
dynamics it would be preferrable to use the fermionic field theories that sup- 
port the static vortex solutions. The simplest fermionic Chern-Simons theory 
that supports the vortex solutions is the one proposed by Li and Bhaduri W, 
described by the Lagrangian, 

C = l e ^ Ffiu A a + i&fW ~ m*H> + eA^f, (34) 

where ip is a two-component spinor field, the Dirac matrices are chosen to be 

7 °=a 3 , 7 1 =ia 1 ,7 2 =^ 2 (35) 
in terms of the Pauli matrices a 1 , and 

f = fyfip, ip = V f 7°. (36) 

We will denote the charge density as j° = ifi'if) = p. 

To find the static vortex solutions in this system, it is simpler to start from 
the field equations, 

j^iidfj, + eA M )V> - mip = , 

ler a F va = -ef . (37) 

The second equation of (|3^) can be decomposed into the Gauss' law constraint, 

B = F 12 = —p, (38) 

K 

and 

E l = F 0l = -«''/. (39) 

K 
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By writing 



the fermion field equation, the first of (|37|), reduces to the two coupled equations 
for ip + (x, t) and ip- [x, t) : 

(id - m)ip + (x,t) = (Di - iD 2 )tp-(x,t) 
-{id Q + m)ip-{x,t) = (D 1 + iD 2 )^+(x,t), (41) 

where Di — di — ieA 1 . 

We now seek the stationary solution of the form, 

^ = ^+(2) ( I ) e~ iEft , (42) 
where Ef is a constant. Then the fermion field equation reduces to 

(E f -m)1>+{g)=0, (43) 
which determines the constant Ef = m, and the self-dual equation, 

{D 1 +iD 2 )^ + {x)=Q. (44) 

Note that this is the same type of self-dual equation as that of the Jackiw-Pi 
model. 

If we take , 

*P = *P-(x) ( J \e~ iE f\ (45) 

instead, the fermion field equation reduces to 

(E f +m)ip-(x)=0 (46) 

which determines Ef = —m and the self-dual equation, 

(D 1 - iD 2 )tp-(x) = 0. (47) 

To find the static solutions we choose the gauge, A = 0, and take A to 
be static. Then the both sides of Eq.(|39|) consistently vanish, and the spatial 
components of matter current vanish for spinor fields (pl2| ) and (|4^), since 7*'s are 
off-diagonal. Thus for the upper component spinor field (|4^) , the field equations 
reduce to the self-dual equation (^4]) and the Gauss' law constraint (|38"|). If we 
write the spinor field as 

i; + (x) = pie i ", (48) 



Eq.(44) reduces to 



eA + = - l -d+ In p + d+tu, (49) 
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which determines the gauge field in terms of the matter fields 

1 



A = Vw V x In p. 

2 1 

From Eq. (f49|) together with the Gauss' law constraint 
matter charge density satisfies the Liouville equation: 



V 2 lnp 



2e 2 



(50) 

one finds that the 
(51) 

k > is required in 



Thus if we take the upper-component for spinor field (|4S 
order to have non-singular positive charge density p. 

If we take the lower-component for the matter field, Eq(45), on the other 
hand, the field equations reduce to 



V 2 Inp 



2e 2 



-(>■ 



(52) 



where k must be negative for the regularity of the charge density p. 

Since the matter charge density p satisfies the same Liouville equation ( |5l| ) 
and ([32]) as that of Jackiw-Pi model, we find the same structure of the static 
vortex solutions except that the energy of the solution is now given by 



E = ±m I d xp± 



(53) 



where p± = tj)± {x)ip±{x) for the upper and lower component matter fields, 
respectively. We thus find that the fermion Chern-Simons theory ( |34| ) supports 
the finite energy static vortex solutions with quantized magnetic flux and charge. 

In general field theories with Chern-Simons term, parity is known to be 
violated by the Chern-Simons term. One can restore the parity invariance by 
introducing appropriate parity partners for each field in the theory. A parity 
invariant fermion Chern-Simons theory is proposed by Hagen JlQ|| . This parity 
invariant fermion Chern-Simons theory also supports the static vortex solutions 
as solutions of Liouville equation . 



IV. Vortex Solutions of a Fermion Maxwell-Chern- 
Simons Theory 

To find a field theoretic model which supports the vortex solutions with phys- 
ically more interesting properties, we introduce a Maxwell term to the gauge 
field part of the Lagrangian and a couple of new interaction terms Jl2[ : 



- ~e a ^F a0 A 1 + iWdrf- m# + + IG^) + \g{^f , 

(54) 
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where Fp V — dpA v — dvA^, J' 1 = ip^f^ip and the new current G M is defined by 

G" = t^PduJp. (55) 
The new gauge coupling term can be written as 

= d^A^PJp + surface term 

= F^e^fJp + surface term, (56) 

which is the magnetic moment coupling in 3-dimcnsional space-time. This is the 
reason why this new coupling is called an anomalous magnetic moment coupling. 

If the gauge field part of the Lagrangian ((HJ) is eliminated, then the theory 
reduces to the 3-dimensional Gross-Neveu model. If we let I — ► 0, g — ► and 
e 2 — * oo with fixed, the theory becomes the fermion Chern-Simons theory 
discussed in the last section. 

Field equations of the new theory are 

d v F» v + ^ vp F vp = -e(J" + IG») 
7"(i0p + eAf^ip - rr# + g{^)4> - el&^{d v Ap) lp ^ = 0, (57) 

the first of which is the field equation for a massive gauge field in 3-dimensions. 

To find static solutions we choose the gauge, ^4° = 0, and take A to be static. 
As in the last section we write the spinor field in component form, 

m 

where Ef is a constant. Then the spinor field equation reduces to the coupled 
equations, 

[E f -m + g(\i/j+\ 2 - |V- 1 2 ) - ele^djAi] ^+ = 
[ - E f - m + g(\iP+\ 2 - \iP-\ 2 ) + eleVfyAi}^ = (59) 

where D± = ± iD 2 , and Dj = di — ieA 1 . 
If we take the upper-component spinor field, 

tf = lM*)( J )e~ iEft , (60) 

then Eq.(|5^) reduces to 

[Ef -m + gp+ - ele io djAi]tl; + = 

D+iP+ = 0, (61) 

where p+ = J+° = \ip+\ 2 is the charge density of the matter field. For this 
choice of the spinor field, we find, 

J+ l = +7 V+ = 

G° = (62) 
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since 7*'s are off-diagonal matrices. Then the gauge field equation reduces to 

e^Fn = --p+ 

K 

djF ij = elG 1 = ele ij djp + . (63) 

For these two equations to be consistent, the coupling constants I and k must 
be related by 

1 = 4 (64) 

This shows that, for the fermion Maxwell-Chern-Simons theory (|54|) to have 
consistent static field equations, one need to introduce the anomalous magnetic 
moment coupling term. 

To find static solutions, therefore, one has to solve Eqs. ( |6l|) and (|63|), which 
now reduce to 

(E f -m + gp + + -e ij d j A 4 ]'4) + = 



D+^+ = (65) 



B = -F 12 = -p. 

K 



Note that, for the first equation of Eq.(|66[) to be consistent with the other two 
equations, one need to require, 



E f = m 

g = ~. (66) 



In other words, in order for the theory ( p4| ) to have consistent static self-dual 
solutions, the quartic coupling constant of the matter field must be determined 
by the Chern-Simons coupling by ( p6[) as in the case of the Jackiw-Pi model. 
To solve the second equation of (|66|), we write 

i>+=p+ 1 *e i ". (67) 
Then this equation determines the gauge field in terms of the matter field: 

eA + = - l -p+ 1 d+p + + d+uj 

eA_= l -p- 1 d_p + + d_u J , (68) 

from which one finds, 

d_A+ - d+A_ = 2iB = --d+d_ lnp+. (69) 

e 
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By substituting this equation to the Gauss' law constraint, the third equation 
of (ptf), one finds the Liouville equation for the matter charge density: 

V 2 lnp + = -— p + , (70) 

K 

where K must be positive for the regularity of the matter density. This shows 
that the fermion Maxwell-Chern-Simons theory ( fiij ) also supports the static 
vortex solutions as solutions of Liouville equation. 
If we choose the lower component spinor field, 

^ = ^-{x)(\\e- iE f\ (71) 

the field equations reduce to 

(E f +m + gp- + -e^djAi)^ = 

D-tP- = (72) 
B=-p-, 

K 

where p- — ip> _%/)-. For these equations to be consistent, we also need to 
require, 



Ef = —m 

9 = --,- (73) 



By writing ip- — p_ 2 e lu) , we again obtain the Liouville equation for the matter 
density 

V 2 lnp_ = -— p_, (74) 

K 

where k < is required for the regularity of the density p_ . 

We now consider the energy of the static vortex solutions. The Hamiltonian 
density is given by 

H = ±m P± + ]-F 12 2 - \g P± 2 + — p± 2 . (75) 



Due to the consistency requirements (|34|), (pq), and (|73|), the last three terms 
of ( |75| ) cancel out, and the energy of the static solutions reduces to 

E = ±m J d 2 p± = ^mN±, (76) 

where 2N± denotes the flux number of the solutions. The cancellation of the 
quartic terms in the Hamiltonian density is also reminiscent of that in the 
Jackiw-Pi model. 
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One may wonder why one studies such a complicated model with an unusual 
interaction term if it gives the same static vortex solutions as those of the Jackiw- 
Pi model or the simple fermion Chern-Simons theory. Although the theory ( |54| ) 
gives the same static solutions as those of the Jackiw-Pi model, the moduli space 
dynamics of this model is quite different from the Jackiw-Pi model because of the 
Maxwell term in the Lagrangian. This may give more interesting, and hopefully 
more realistic vortex dynamics. 

The theory we have described is a U(l) gauge theory. This theory can be 
generalized to SU (N) gauge theories as has been done by Jackiw and Pi H for 
the Jackiw-Pi model. There exist many possibilities to formulate corresponding 
SU(N) gauge theories as a generalization of the fermion Maxwell-Chern-Simons 
theory. The simplest example is to take adjoint representation for the matter 
field and use the ansatz that the matter density and the gauge field lie in 
the Cartan subalgebra of the grhoup Q . One then finds the equations satisfied 
by the components of matter density, 

2e r 

S7 2 \np a = T—J2K Q ppf 3 , (77) 

0=1 

where K a p is the Cartan matrix. Eq . (J77| ) is not completely integrable in general, 
but can be integrated numerically For the SU(2) case, Eq.([77|) reduces to 
the Liouville equation, and it corresponds to an embeding of U(l) in SU(2). 



V. Discussions 



We have constructed a fermion Maxwell-Chern-Simons theory that possesses 
completely integrable multi- vortex solutions, and have shown that the solutions 
are the same as those of Jackiw-Pi model. Although the model gives the same 
static solutions as the Jackiw-Pi model, the moduli space dynamics of the static 
solutions will be quite different due to the Maxwell term in the Lagrangian. The 
Maxwell term is quadratic in the time-derivatives of gauge field, i.e., ^AiAi, 
which gives rise to a non-trivial contribution to the moduli space metric, while 
the Chern-Simons term is linear in the time-derivatives of gauge field. This may 
give an interesting vortex dynamics and may be more relevant in understanding 
the vortex phenomena arising in the high temperature superconductors. 

Similar construction can be done for more relatistic, non-relativistic fermion 
gauge theories. Some models have already been constructed and their static 
solutions have been studied M. These models provide simple starting point in 
studying more realistic vortex dynamics in the field theoretic framework. 

The most well-known method of studying vortex dynamics is that of mod- 
uli space dynamics for slowly moving vortices proposed by Manton [ p^[ . This 
method, however, has a difficulty in obtaining the moduli space metric for the 
theories with kinetic terms linear in time-derivatives, such as Chern-Simons the- 
ory, relativistic fermion theories and general non-relativistic field theories. This 
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difficulty can be avoided by using the method proposed by Bak and Lee fll5[ , 
where one integrates out the momentum variables to make the kinetic terms 
quadratic in time-derivatives of field variables. One can then use the Manton's 
method to obtain the correct moduli space metric. 

Another way to study the vortex dynamics is the so-called string formulation 
of vortex dynamics Jl6) . The authors of ref . jl6| note that at the centers of vor- 
tices the field function(scalar field <j> in the Jackiw-Pi model and spinor field ip 
in the fermionic models) vanishes, and treat the points of zeros in 3-dimensional 
space-time (or string of zeros in 4-dimensional space-time) as fundamental ob- 
jects. They then write exact equations of motion for points(or strings) in terms 
of fields that surround them. These equations of motions are rather complicated 
and need approximation for practical applications. 

As mentioned above, it is not a simple matter to use these methods to study 
vortex dynamics in the field theoretic framework. The field theoretic models 
discussed here have static vortex solutions in closed form, and it is relatively 
simple to apply these methods to understand the dynamics of magnetic vortices. 
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